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Проведён анализ сущности алгоритмов вероятностного объединения данных в случае сопровождения одной или нескольких це- 
лей, когда неизвестно, какие измерения порождены целью, а какие появились вследствие переотражений сигнала. Показан ре- 
зультат работы алгоритмов и приведены статистические свойства оценок. 


Введение 

На практике часто возникает задача определе- 
ния параметров движения объектов по локацион- 
ным измерениям. Эту задачу приходится, напри- 
мер, решать при определении орбит ИСЗ, а также 
при управлении движением судов и самолетов дис- 
петчерскими службами морских, речных портов и 
аэропортов. В данной работе рассматриваются ме- 
тоды определения параметров движения объектов 
в статистически неопределенных ситуациях, при 


которых неизвестны соответствия между измере- 
ниями и объектами, возможен пропуск измерений, 
появление ложных измерений. 

Задача сопровождения одиночной или группо- 
вой цели при наличии ложных сигналов характер- 
на для многих типов радиосистем. В условиях, ког- 
да среди принимаемых сигналов могут быть как из- 
мерения от цели, так и ложные тревоги, возникает 
проблема отождествления отсчетов наблюдений, 
которые следует использовать для получения оце- 
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нок координат в измерительных радиосистемах. 
Для решения подобной задачи, применительно к 
сопровождению одиночной цели, широко исполь- 
зуется метод вероятностного объединения данных 

[1] . Впервые о нем упоминается в 1970 г. Затем по- 
явились его многочисленные модификации: сов- 
местное вероятностное объединение данных, ста- 
тистически упорядоченное совместное вероят- 
ностное объединение данных и др. 

1. Методы сопровождения 

В случае сопровождения целей в условиях, ког- 
да существуют ложные цели и помехи, возникает 
проблема - какое из совокупности принятых изме- 
рений нужно использовать для обновления каждой 
траектории. Существует множество методов, соз- 
данных для решения этой задачи. Два простых под- 
хода: «сильнейший сосед» и «ближайший сосед» 

[2] . В методе «сильнейший сосед» используется 
сигнал с наибольшей мощностью из всех измере- 
ний, попавших в строб селекции (строб селекции - 
область вокруг предсказанного измерения, исполь- 
зуемая для выбора измерений, которые будут уча- 
ствовать в обработке), остальные измерения отбра- 
сываются. В методе «ближайший сосед» использу- 
ется измерение, расположенное ближе всего к сво- 
ему предсказанному значению. Эти простые мето- 
ды работают хорошо, когда отношение сигнал/шум 
велико и цель не маневрирует, и теряют свою эф- 
фективность, когда увеличивается норма ложных 
целей или когда цель маневрирует. 

Помимо использования только одного измере- 
ния среди полученных и отбрасывания остальных, 
существует и альтернативный подход, который ис- 
пользует все измерения, попавшие в строб селек- 
ции, но с разными весами (вероятностями). Этот 
метод называется вероятностным объединением 
данных. Данный алгоритм и его многочисленные 
улучшенные версии показали большую эффектив- 
ность в сопровождении цели на фоне помех. 

Проблема объединения данных становится 
сложнее при сопровождении нескольких целей, 
когда траектории «конкурируют» между собой за 
измерения. Здесь в дополнение к тому, что каждая 
траектория может иметь в области строба несколь- 
ко измерений (как и в случае одиночной цели), так 
и само измерение может попадать в стробы сразу от 
нескольких траекторий и тогда идёт «соревнова- 
ние» между траекториями за измерение. Существу- 
ет множество алгоритмов для разрешения данного 
противоречия, например, метод «глобального бли- 
жайшего соседа» (создаёт «жесткое» распределение 
измерений по траекториям) и метод совместного 
вероятностного объединения данных [3]. 

2. Вероятностное объединение данных 

Алгоритм вероятностного объединения данных 
используется для сопровождения одиночной цели 
в присутствии ложных целей и помех. Его суть со- 


стоит в том, что каждое из измерений, попавшее в 
строб селекции, по очереди принимается за истин- 
ное. С использованием соотношений стандартного 
фильтра Калмана [4] формируются оценки, кото- 
рые затем объединяются в одну с весами, соответ- 
ствующими вероятности того, что данное измере- 
ние порождено целью. 

Данный алгоритм строится исходя из следую- 
щих предположений [ 1 ] : 

• сопровождается только одна цель, траектория 
которой уже завязана; 

• динамика л-мерного вектора состояния цели 
х(к) на к- ом временном такте задана матрич- 
ным разностным уравнением: 

х(к + 1) = Р(к)-х(к), 

где Ѵ(к) — матрица перехода; 

• линейная модель наблюдений определена соот- 
ношением 

г(к) = ЩА:) • х(к) + уѵ(к), (1) 

где Н(^) - матрица наблюдений, ѵ/(к) - гауссова 
последовательность белого шума с нулевым сред- 
ним и известной ковариационной матрицей; 

• информация о цели к моменту времени к с уче- 
том всех данных до этого момента Ъ к ~ 1 \ 

р[х(к)/% к - 1 ] = М[х(к);х(к/к-1), Р(к/к-1)], 

где р[х{к)/Ъ^ ѵ \ - условная вероятность; 

Щх(к)',х(к/к-\),Р(к/к-1)] - гауссова плотность ве- 
роятности с аргументом х(к), средним значением 
х(к/к— 1) и ковариационной матрицей ошибки 
оценки Р(к/к-\); 

• на каждом такте в фильтре сопровождения фор- 
мируется строб по сигналу (1). При вероятности 
обнаружения Р в отметка от цели может оказать- 
ся внутри строба с вероятностью Р 0 \ 

• в строб попадает не более одной отметки, по- 
рождённой целью. 

Исходные данные: задаются Р в , Р 0 , х(0/0), ДО/О). 
Получив оценку состояния на к- ом такте и рас- 
считав ковариационную матрицу ошибки оценки, 
обращаются к уравнениям предсказания на сле- 
дующий такт: 

х(к + \/к) = Р(к)х(к/к), 

Р {к + 1/к) = Р(к)Р(к/к)Р т (к). 

Затем алгоритм начинает работу заново, т. е. яв- 
ляется рекурсивным. 

Недостатком данного алгоритма является то, 
что количество вычислений экспоненциально ра- 
стёт с числом измерений. Но при приемлемом 
уровне интенсивности помех он является достаточ- 
но эффективным в случае сопровождения одиноч- 
ной цели. 

В ходе выполнения данной работы были иссле- 
дованы статистические свойства оценок, получен- 
ные с помощью алгоритма вероятностного объеди- 
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нения данных. Для этого была разработана про- 
грамма в пакете МаІІіСАО и методом прямого ве- 
роятностного моделирования (метод Монте-Кар- 
ло) получены оценки дисперсий ошибок. 

Двумерный вектор состояния динамической 
системы имеет вид: 



а-і+Ь ' 
а ) 


( 2 ) 


где а - скорость изменения дальности до цели; Ъ - 
постоянная составляющая (начальное значение 
дальности). Для измерения доступна только первая 
координата х г (() вектора состояния, т. е. мгновен- 
ное значение дальности до цели. Матрица перехода 
имеет вид 


Г(к) 



где А і - шаг дискретизации. 


На рис. 1 изображена динамика изменения ве- 
личины строба, оценок состояния и наблюдаемых 
отметок для случая, когда среднее число ложных 
отметок в стробе в начальный момент времени рав- 
но 3,7, отношение сигнал/шум 

<? = <[Е х і 2 (к) 

V *= і 

где N=20 - количество отсчетов на временной оси, 
а ш - среднеквадратическое отклонение (СКО) шума. 

Исследования СКО ошибки сопровождения про- 
изводились усреднением по 1500 независимым реа- 
лизациям оценок. На рис. 2 приведена зависимость 
величины <7„ (СКО оценки, нормированное к своему 
априорному значению) от времени для различных 
средних значений числа ложных отметок в стробе и 
для значений отношения сигнал/шум #=100 и #=50. 

Расчеты показывают, что с ростом среднего чи- 
сла ложных отсчетов увеличивается СКО ошибки. 




Рис. 1. Динамика изменения величины строба, оценок состояния и наблюдаемых отметок (одиночная реализация) 



к к 

Рис. 2. Нормированная СКО ошибки сопровождения 
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В исследуемом диапазоне среднего числа ложных 
отметок в стробе с=1,5...6,0 фильтр сопровождения 
на основе алгоритма вероятностного объединения 
данных обеспечивает снижение СКО оценок по 
сравнению с априорным значением в 10. ..2, 5 раз 
при <7=100 и в 3,5. ..2,5 раз при <7=50. 

При работе алгоритма возможны случаи, когда 
происходит его срыв. Если интенсивность помех 
велика, то возможна ситуация, при которой лож- 
ное измерение оказывается ближе к предсказанно- 
му, т. е. оно имеет больший вес при формировании 
оценки. В результате получается неверная оценка, 
вокруг которой затем выставляется новый строб, и 
понятно, что истинное измерение еще с большей 
вероятностью окажется дальше от своего предска- 
занного значения. Этот процесс имеет лавинооб- 
разный характер. Оценка на выходе фильтра все 
дальше уходит от истинного измерения, которое 
все слабее (с малым весом) влияет на формирова- 
ние очередной оценки. В итоге она формируется в 
основном по ложным отсчетам, и наступает срыв 
слежения. Негативно влияет на работу фильтра и 
уменьшение отношения сигнал/шум. В частности, 
при <7=50 и с= 3,7 оценка вероятности срыва равна 
р=0,014; при <7=100 она равна р=0,012, а при <7=50 и 
с=6: 77=0,023. Эффективность работы фильтра су- 
щественно зависит от отношения сигнал/шум и 
интенсивности помех. 

3. Совместное вероятностное объединение данных 

Вероятно, наиболее сложной проблемой в тео- 
рии обработки радиолокационной информации 
является сопровождение нескольких целей в при- 
сутствии помех. Эта проблема возникает, если цели 
находятся на небольших расстояниях друг от друга, 
что, в частности, имеет место при пересечении их 
траекторий или при полёте группой. В этих случаях 
стробы селекции накладываются друг на друга, и 
получаемые отметки могут быть отнесены к нес- 
кольким траекториям. Фильтр на основе вероят- 
ностного объединения данных может использо- 
ваться в данной ситуации, но его возможности 
очень ограничены, т. к. модель помех (равномерное 
распределение в стробе селекции) больше не отве- 
чает действительности. 

Поскольку истинное измерение обнаруживает- 
ся, как правило, с вероятностью Р в Р е > 0,7, то его 
наблюдение можно считать относительно устойчи- 
вым на протяжении нескольких тактов наблюде- 
ния. Ложные измерения, с другой стороны, дол- 
жны быть независимыми и с пренебрежительно 
малой вероятностью способны образовывать 
устойчивые последовательности от такта к такту 
(если интенсивность помех высока). Таким обра- 
зом, можно допустить, что в области выставленно- 
го на каждом ^-ом шаге строба, они распределены 
равномерно. Если области стробов от двух и более 
целей пересекаются, то это может означать появле- 
ние устойчивых общих для этих целей источников 
помех, что, конечно, ухудшает характеристики 


фильтра, основанного на алгоритме вероятностно- 
го объединения данных. В такой ситуации одно и 
то же измерение может принадлежать различным 
объектам или быть ложным. 


Совместное вероятностное объединение дан- 
ных [3] решает задачу получения оценки на основе 
вычисления вероятности объединения совместно 
для всех траекторий целей, чьи области стробов пе- 
ресеклись. При этом рассматриваются все возмож- 
ные привязки «измерение-траектория» с расчётом 
вероятностей этих привязок. Итоговая оценка 
представляет собой взвешенную сумму оценок по 
каждой из привязок. 

Алгоритм опирается на два важных предполо- 
жения: ни одна отметка не может быть порождена 
более чем одной целью, и за один обзор никакая 
цель не должна создавать более одной отметки. 
Данный алгоритм математически отличается от ал- 
горитма вероятностного объединения расчётом ве- 
совых коэффициентов (3/(к ), равных вероятности 
того, что у-ое измерение принадлежит цели I. Он 
эффективно работает в случае сопровождения нес- 
кольких целей, но обладает таким же недостатком, 
как и алгоритм вероятностного объединения дан- 
ных - объём вычислений растёт экспоненциально 
с числом целей и числом ложных отметок. 


Существует алгоритм линейного вероятностно- 
го объединения данных. Он тоже используется для 
сопровождения групповой цели. Но в отличие от 
алгоритма совместного вероятностного объедине- 
ния данных, где число операций растёт экспонен- 
циально с числом целей, в этом алгоритме данная 
зависимость линейна и при этом наблюдается 
лишь незначительное ухудшение характеристик 
сопровождения. Объём вычислений уменьшен за 
счет того, что при рассмотрении одной цели, дру- 
гие траектории обрабатываются в качестве допол- 
нительных источников шума. Это позволяет со- 
провождать групповую цель, используя вычисли- 
тельные ресурсы алгоритма сопровождения оди- 
ночной цели. Применение данного алгоритма сов- 
местно с интерактивной мультимоделью (набор 
фильтров, выполняющих алгоритм объединения 
данных, на каждую из моделей движения цели) по- 
зволяет сопровождать маневрирующие цели [5]. 

В случае большой плотности помех при сопро- 
вождении нескольких целей, когда алгоритм сов- 
местного вероятностного объединения данных уже 
не справляется, эффективно использовать алго- 
ритм статистически упорядоченного совместного 
вероятностного объединения данных. Он состоит в 
следующем: весовые коэффициенты в алгоритме 
совместного вероятностного объединения данных 
можно записать в виде матрицы 


в = {р)} 


Рі РІ - Рі 
А 2 Рі - РІ 

_р: рі - р:_ 
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Эта матрица упорядочивается: наибольшее зна- 
чение выставляется в первом столбце, а наимень- 
шее - в последнем: 


в ш = Щ} = 


Р 

Р 


1 

Ш,тах 

2 

1.57, тах 


Р 

Р 


1 

2 нс/, шах 
2 

2 псі, шах 


.1 

тік,т ах 

.2 

тік,т ах 


Р 


Т 

І5/,тах 


Р 


Т 

2т/,тах 


Р, 


Т 

тік, тах 


Матрица (3) составляется, чтобы определить 
доминантное измерение. Доминантное измерение 
первой цели - то измерение, которое имеет макси- 
мальное значение весового коэффициента. Доми- 
нантное измерение цели 1 выбирается таким обра- 
зом, что его весовой коэффициент максимален. Но 
не является доминантным для целей с номерами от 
1 по Г-1. Формируется новая матрица, в которой 
весовые коэффициенты, не соответствующие до- 
минантным измерениям для цели і, берутся с мас- 
штабным коэффициентом меньше единицы. Дан- 
ный алгоритм обеспечивает улучшение точностных 
характеристик сопровождения. 


Работа данных алгоритмов была смоделирована 
в среде МаіЬаЪ. Дальность до цели изменяется по 
линейному закону (2). На рис. 3 изображена общая 
картина отметок от цели и результат их обработки. 

В ситуации, когда дальность до разных целей 
отличается на небольшую величину, показательно 
проиллюстрировать работу алгоритмов на рис. 4. 

На рис. 4 жирной линией отмечена оценка из- 
менения дальности, а тонкой пунктирной - истин- 
ное изменение дальности. Видно, что при исполь- 
зовании алгоритма совместного вероятностного 
объединения данных происходит смешение оце- 
нок, т. е. две цели принимаются за одну, вторая 
цель срывается с сопровождения. В этой же ситуа- 
ции алгоритм статистически упорядоченного сов- 
местного вероятностного объединения данных по- 
зволяет сопровождать две цели раздельно. 

Заключение 

Рассмотрены основные методы сопровождения 
целей на фоне помех. Задача становится достаточно 



Рис. 3. Общая картина отметок: а) полученных в течение всего времени обработки и б) формирование по этим отметкам оцен- 
ки дальности до трёх целей 



Время 



Время 

б 


Рис. 4. Оценка дальности до двух целей с помощью алгоритма совместного вероятностного объединения данных: а) стандарт- 
ного и б) статистически упорядоченного 


71 




Известия Томского политехнического университета. 2006. Т. 309. № 8 


сложной, когда число объектов неизвестно. В слу- 
чае сопровождения цели при наличии ложных сиг- 
налов следует вычислять оценку не по одному из 
принятых отсчётов, отвергая все остальные, а на 
каждом такте уточнять оценку интересующего па- 
раметра по всем измерениям (с весами, равными 
апостериорным вероятностям принадлежности из- 
мерений сопровождаемому объекту). 

Методом прямого статистического моделирова- 
ния выполнен анализ работы фильтра вероятност- 
ного объединения данных. В исследуемом диапазо- 
не среднего числа ложных отметок в стробе фильтр 
сопровождения обеспечивает снижение СКО 
ошибки оценки по сравнению с априорным значе- 
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нием до 10 раз. Эффективность работы фильтра за- 
метно зависит от отношения сигнал/шум и интен- 
сивности помех. 

Алгоритм вероятностного объединения дан- 
ных, работающий с одиночной целью, является ос- 
новой для формирования различных его модифи- 
каций. Например, в случае групповой цели, когда 
стробы, выставленные для отдельных объектов, пе- 
ресекаются, необходимо использовать совместное 
вероятностное объединение данных. Недостатком 
алгоритма вероятностного объединения является 
экспоненциальный рост вычислений с увеличени- 
ем интенсивности ложных отметок и количества 
сопровождаемых целей. 


4. Сейдж Э. Теория оценивания и её применение в связи и упра- 
влении. - М.: Связь, 1976. - 320 с. 

5. Мікісі Б., 8иѵогоѵа 8., Сііаііа 8. Миііі Іаг§еІ Ігаскіп§ оГ §гоипб 
Іаг^еП іп сІиМегмШ ЬМШБА-ММ //ТЬе Т" Іпіегп. СопГ. оп Іп- 
Гогтайоп Ризіоп. - 8іоскйо1т, 2004. - Р. 1 104-11 10. 


ѴПКЮ1 ЗДК РМЗЗ 




